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Abstract 

In this paper we examine the data that was collected at Haley Station in Antarctica 
on June 22, 1987. Using a test devised by Dewan ^ we interpret the flow as one which 
represents two-dimensional turbulence. We also construct a model to interpret the 
spectrum of this data which is almost independent of the wave number for a range of 
frequencies. 
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1 Introduction 



Two dimensional turbulence has been the subject of intense theoretical research [||, |[ and 
simulation experiments The reason for this interest stems from the fundamental differ- 
ences between 3-d isotropic and 2-d turbulence. To begin with, vortex stretching is absent 
in 2-d as a direct consequence of Navier-Stokes equations. Furthermore in 3-d the energy 
cascade is from the large eddies to small one but this process reverses itself in 2-d and leads 
to the formation of large scale coherent eddies. Another difference between two and three 
dimensional turbulence exists in the inertial range of the spectrum. Kraichnan showed Q 
that in 2-d in addition to Kolmogorov inertial range there is (due to ensotrophy conservation 
in zero viscosity) another scaling law in the form 

where t] is ensotrophy dissipation rate. 

While many simulations P, ^ confirm these theoretical predictions the actual observation 
and detection of 2-d turbulence as a natural phenomena remains (as far as we know) an open 
questions. 

One of the objectives of this paper is to weigh in the pros and cons for 2-d turbulence 
in the Antarctic data that was obtained by the British observation post as Haley Station in 
Antarctica on June 22, 1987 (for further description of this data see 0, H). The importance of 
these measurements stem from the fact that the flow field u = {u, v, w) and the temperatures 
were measured simultaneously at three different heights viz. 5m, 16m and 32m. These 
simultaneous readings enable us to apply a test devised by E. Dewan ||^ for the detection of 
2-d turbulence. According to this test 2-d turbulence is characterized by small values for the 
coherence |^ between the time series which represent the various meteorological variables 
at different heights. 

From another point of view the Antarctic data represent a stably stratified medium. 
(According to mission records the temperature gradient with height can reach up to IK/m). 
Under these circumstances Bolgiano [0, [111 and others p speculated about the existence of 



"buoyancy range turbulence" (BRT) which should lead to a flattening of the spectra in parts 
of the inertial range. In this paper we shall estimate the power spectrum for the data using 
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the usual Fourier transform and by the method of maximum entropy (briefly the reason 
for this duplicatcy is due to the existence of "discontinuities" in the data). Both of these 
estimates show a spectral range in which the spectrum is almost flat and thus support the 
theoretical arguments that were advanced for the existence of BRT. 

The plan of the paper is as follows: In section 2 we describe the method used to filter 
out the mean flow and waves from the data and the tests that were applied to verify that 
the residuals actually represent turbulence. In section 3 we apply the coherence test for 
2-d turbulence and discuss its consequences. In section 4 we present a model for the power 
spectrum of the data and its implications. We end up in section 5 with some conclusions. 

2 Data Detrending 

The statistical approach to turbulence splits the flow variables u, T (where T is the temper- 
ature) into a sum 



where u, T represent the mean (large scale) flow, u\ T' represent waves and m, t "turbulent 
residuals" [[g 

To effect such a decomposition in our data we used the Karahunan-Loeve (K-L) decom- 
position algorithm (or PCA) which was used by many researchers (for a review see [0). 
Here we shall give only a brief overview of this algorithm within our context. 

Let be given a time series X (of length A^) of some geophysical variable. We first deter- 
mine a time delay A for which the points in the series are decorrelated. Using A we create 
n copies of the original series 



(To create these one uses either periodicity or choose to consider shorter time-series). Then 
one computes the auto-covariance matrix R = (Rij) 



u = u + u' + u, T = T + T' + t 



X{k), X{d + A),... ,X{k + {n -1)A). 



N 




(2.1) 



k=l 
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Let Ao > Ai, . . . , > A„_i be the eigenvalues of R with their corresponding eigenvectors 



^o>--- ,0n-i)> i = 0,... ,n- 1. 



The original time series T can be reconstructed then as 

n-1 

X{j) = Y,<^k{M (2.2) 

k=0 

where 

^ n— 1 

a.(j) = -E^(.^' + ^^)<^'- (2.3) 



i=0 



The essence of the K-L decomposition is based on the recognition that if a large spectral 
gap exists after the first mi eigenvalues of R then one can reconstruct the mean flow (or 
the large component ( of the data by using only the first mi eigenfunctions in ( p.2|) . A 
recent refinement of this procedure due to Ghil et al |jl3| is that the data corresponding to 
eigenvalues between mi + 1 and up to the point m2 where they start to form a "continuum" 
represent waves. The location of m2 can be ascertained further by applying the tests devised 
by Axford fl^ and Dewan [^] (see below). 



Thus the original data can be decomposed into mean flow, waves and residuals (i.e. data 
corresponding to eigenvalues m2 + 1, . . . , n — 1 which we wish to interpret at least partly as 
turbulent residuals). 

For the data under consideration we carried out this decomposition using a delay A of 
1024 points (approximately 51 sec.) for all the geophysical variables. In table 1 we present 
the values of mi, m2 that were used in this decomposition for the flow variables at different 
heights. (In all cases n = 64). 

The residuals of the time series which are reconstructed as 

n-1 

X'~{j)= «^(^>o (2.4) 

k=m2+l 

contain (obviously) the measurement errors in the data. However to ascertain that they 
should be interpreted primarily as representing turbulence we utilize the tests devised by 
Axford and Dewan 0. According to these tests turbulence data (at the same location) 



is characterized by low coherence between u,v,w and a phase close to zero or vr between w 
and t. (A phase close to 7r/2 is characteristic of waves). Figs. 1,2,3 show samples of the 
coherence between the residuals of u, v, w at different heights. They demonstrate that for 
most frequencies the coherence is less than 0.1. Fig. 4 gives a scatter plot of the phase 
between w and t at height 5m. This figure is less definitive as there are still quite a few 

— , ]■ However out of the 200 points in this plot 125 are in 

the "turbulence sector". 

These tests show that to a large extent the residuals that were obtained from the K-L 
decomposition represent actual turbulence. 



3 Tests for 2-d turbulence 



In today literature [|Ig] a spectral slope of —3 in part of the inertial range is considered to 
be a strong indicator for 2-d turbulence. However as noted already by Lily |^ "geophysical 
consideration" might modify this slope. Since the spectral plots for the flow under consider- 
ation (for sample see figs. 8,9,10) do not exhibit this dependency (except for w at 16m in the 
low frequencies) we must resort to other tests to bolster the claim that the flow described 
by this data corresponds to 2-d turbulence. 

To this end we utilize a test devised by Dewan 0. According to this test inviscid two 
dimensional turbulence is characterized by the fact that the temporal statistical coherency 
^ between the time series representing the flow variables at different altitudes is zero. With 



viscosity taken into account some vertical separation of the order of (10m for air) is needed 
for the coherency to become small. (Strong coherency with values close to one indicates a 
strong linear relationship between the two time series pO[| ). 



Some typical plots for the coherency in the data is presented in figs. (5,6,7). In these 
plots the coherency for w between the different heights is plotted for different wave numbers. 
We observe that for most sampled frequencies the coherency is well below 0.1 and according 
to Dewan "these values constitute evidence for 2-d turbulence and against other types of 
fluctuations" . 
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4 The spectrum 



Two dimensional flow of incompressible and inviscid fluid conserve both the energy E and 
the enstrophy Q. Fir viscous fluid these quantities decay according to 



dE 



-e = — = -2uQ, -e^ = — = -u \ Vuj 
ot ot 



(4.1) 



The energy spectrum is determined therefore by both parameters e, which leads to the 
definition of a length scale 

1/2 



L.. 



(4.2) 



From dimensional considerations one concludes then that |]T6[ the energy spectrum in the 
inertial range must have the form 



(4.3) 



where / is a function of the dimensionless variable fcL^. If at one end of the inertial range only 
e is essential (and the effect of is negligible) then / = constant and the energy spectrum 
obey Kolmogorov 5/3 power law. If on the other end of this range e is not essential then / 
must have the form 



and consequently 



E{k) = Ce'fk-' 



(4.4) 



(4.5) 



(where C is a constant). 



For stratified medium Obukov |17| introduced the temperature inhomogeneity dissipation 



rate 



POD 

er = 2x / k'^ET{k)dk 
Jo 



(4.6) 



where Et is the temperature spectra and x is the heat conductivity of the medium. He 
further postulated that the turbulent component of T is dependent on this parameter. 
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For the (stratified) Antarctic medium we would like to enlarge the domain of this pos- 
tulate to include the velocity components of the flow. This enables us to introduce the 
buoyancy (length) scale |l^ 



Lb = {ag)-'/h'/%'/' (4.7) 

where {ag) is the buoyancy parameter. The existence of this second length scale for stratified 
two dimensional flow lead us to replace ([4.3[) by 



E{k) = f{kL^, kLB)e'/'k-'/^ (4.8) 

However since stratification and enstrophy conservation are independent of each other we 
infer that / must have the form 

/ ^ {kL^YikLBY. (4.9) 

It follows then that the spectral dependence on k is given by 

E{k) ~ (4.10) 

We conclude therefore that various combinations of r, s are possible and this will lead to 
different spectral dependencies on k. 
Thus if 

E{k) ~ k-'^ 

and the dissipation e is negligible we must have then 

r 5 2 
r + s = 5/3-g, 2^4*^3"'^ 

which yields 

33 - 15g 15q - 18 

From the spectral plots for the data under consideration we see that (approximately) 

E{k) ~ 

for a large segment of the inertial range which is characteristic of the "buoyancy range 



turbulence" as predicted by Bolgiano |T0|, ITT . 



It is interesting to note in this context that Kriachnen []I9[ already observed that the 
"energy spectrum of the flow depends on the details of the nonlinear interaction embodied 
in the equations that govern the flow and can not be deduced solely from the symmetries, 
invariances and dimensionality of the equations" . 

Finally we would like to observe that the data under consideration contains some dis- 
continuities. These can change completely the asymptotic behavior of the spectrum. To 
demonstrate this assume that the data is described by 



D{x) = CH{x - xo) + g{x) 



(4.11) 



where g{x) is a smooth function whose Fourier transform (FT) decays exponentially and 
H{x) is the Heaviside function 



Hix) 



1 x > 
X < 0. 



Differentiating ( [4.11 ) we have 



D\x) = C6{x - Xo) + g'ix) 



(4.12) 



and the FT of (|l|) is 



D'{k) = C+~g'ik) 



(4.13) 



The FT of D is obtained then by dividing ( [4. 131) by k which shows clearly that the asymptotic 
behavior of D{k) is proportional to k~^. 

We conclude then that a proper filter for the removal of these discontinuities from the 
data is needed in order to obtain the true spectrum of the turbulent residuals. Such a filtering 
algorithm is given by the K — L decomposition which was described in Sec. 2. 



5 Conclusion 

Using the coherency test advanced by Dewan we are able to characterize the flow under 
consideration as one that has the characteristics of 2-d turbulence. One stumbling block for 
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this interpretation is the absence of —3 slope in part of the inertial range. To explain this we 
introduced a model that takes into account the stratification of this flow. This model shows 
that when buoyancy effects are taken into account different slopes of E{k) are possible. Thus 
we believe that we introduced evidence for the interpretation of this spectra as one belonging 
to BTR. 
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m2 


J- f 

u at 5m 
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42 


V at 5m 
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w at 5m 
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i at 5m 
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u at 16m 
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V at Ibm 
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w at Ibm 
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1 at Ibm 
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u at dim 
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V at 32m 
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w at 32m 
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T at 32m 


2 
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Table 1 
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